j 

i 


Reproduced  by  Ihi 

CLEARINGHOUSE 
lor  federal  Scientific  &  Technical 
Information  Springfield  Ve.  2x151 


AVW?/) 


FINITE  PLATE  BENDING  ELEMENTS  IN  POLAR  CO-ORDINATES 


by 

MERVYN  D.  OLSON,  GARRY  M.  LINDBERG,  HELEN  A.  TULLOCH 


A.  H.  Hall,  Head 

Structures  and  Materials  Section 


F.  R.  Thurston 
Director 


SUMMARY 


Two  finite  elements  for  plate  bending  in 
polar  co-ordinates  are  derived.  One  element  is  in  the 
form  of  a  circular  sector  and  the  other  is  in  the  form 
of  an  annular  sector.  The  former  element  has  nine 
degrees  of  freedom,  while  the  latter  has  twelve,  the 
degrees  of  freedom  in  each  case  being  a  displacement 
and  two  slopes  at  each  element  corner. 

The  method  is  checked  with  several  numer¬ 
ical  examples.  The  static  deflections  of  an  annular 
plate  and  a  complete  circular  plate,  bcth  loaded  by 
single  concentrated  forces,  are  analyzed  and  the  re¬ 
sults  compared  with  exact  solutions.  Calculations  of 
the  free  vibrations  of  complete  circular  plates  are 
also  carried  out  and  compared  with  exact  solutions. 
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FINITE  PLATE  BENDING  ELEMENTS  IN  POLAR  CO-ORDINATES 


i.  0  INTRODUCTION 

The  finite  element  method  of  structural  analysis  is  now  firmly  established  as 
a  powerful  technique  for  handling  difficult  problems  in  solid  mechanics.  Finite  elements 
have  been  developed  in  many  shapes  and  for  varied  applications,  the  most  numerous 
being  the  rectangular  and  triangular  plate  elements  for  both  bending  and  plane  stress 
configurations.  A  good  review  of  the  state  of  the  art  can  be  found  in  the  text  by 
Zienkiewicz  (Ref.  1),  but  even  this  book  is  becoming  outdated  by  the  rapidly  advancing 
field. 


The  simplest  element  shapes  for  plate  problems  are  obviously  a  triangle  and 
rectangle  with  three  and  four  nodes,  respectively.  Rectangular  elements  are  somewhat 
limited  in  their  boundary  shape  applicability,  whereas  triangular  elements  may  be  used 
to  represent  almost  any  shape  of  boundary.  However,  for  problems  with  curved  bound¬ 
aries,  this  use  of  triangular  elements  means  that  the  curved  boundary  is  being  approxi¬ 
mated  by  a  series  of  straight  line  segments.  Further,  with  some  of  the  more  sophisti¬ 
cated  elements  being  developed  today,  the  error  introduced  by  this  approximation  may 
well  be  the  limiting  factor  in  some  solutions. 


Hence,  there  appears  to  be  a  need  for  elements  with  curved  boundaries.  A 
start  in  this  direction  has  been  made  by  Ergatoudis  etal.  in  Reference  2,  wherein  quad¬ 
rilateral  plane  stress  elements  with  various  curved  sides  are  developed.  It  would  be 
useful  to  have  this  approach  extended  to  bending  applications  as  well.  In  particular, 
there  is  a  large  class  of  problems  involving  circular  arc  boundaries  that  are  difficult 
to  handle  by  ordinary  analytical  means.  Problems  commonly  encountered  in  jet  engines, 
such  as  the  vibrations  of  turbine  discs  and  turning  vanes,  are  but  two  examples. 


In  the  present  work,  a  method  for  handling  such  problems  is  considered.  Two 
plate  bending  finite  elements  are  developed  in  polar  co-ordinates,  which  may  be  joined 
together  in  various  combinations  with  rectangular  elements  to  fit  a  large  number  of 
boundary  shapes.  The  first  element  has  three  nodal  corners  and  is  in  the  shape  of  a 
segment  of  a  circle.  The  second  element  has  four  nodal  corners  and  is  in  the  shape  of 
a  segment  of  an  annulus.  The  transverse  displacement  and  two  slopes  are  used  as  the 
degrees  of  freedom  at  each  element  nodal  corner.  Displacement  functions  are  derived 
for  each  element,  and  the  stiffness  and  mass  matrices  are  obtained  by  substituting  these 
displacements  into  the  strain  and  kinetic  energy  integrals  from  plate  theory. 


These  elements  are  fully  compatible  with  the  well-known  twelve  degree  of 
freedom  nonconforming  rectangular  plate  elements  (Ref.  1).  That  is,  either  of  the 
polar  co-ordinate  elements  developed  herein  may  be  joined  to  one  of  these  rectangular 
elements  by  equating  the  displacement  and  edgewise  slope  degrees  of  freedom  at  the 
ends  of  a  common  boundary.  The  transverse  displacement  will  then  be  continuous 
across  such  a  junction,  but  the  normal  slope  will,  in  general,  not  be  continuous  there. 
This  is  consistent  with  the  fact  that  these  rectangular  elements  are  net  conforming. 
Hence,  these  polar  co-ordinate  elements,  when  combined  with  the  rectangular  elements, 
provide  a  powerful  method  for  handling  a  variety  of  plate  bending  problems  encountered 
in  practical  applications. 
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In  the  present  work,  numerical  calculations  are  carried  out  only  with  combi¬ 
nations  of  the  polar  co-ordinate  elements.  The  examples  are  chosen  primarily  as  tests 
of  the  elements.  The  static  deflection  of  an  annular  plate  and  a  complete  circular  plate, 
both  loaded  by  single  concentrated  forces,  are  analyzed,  and  the  results  compared  with 
exact  solutions.  Finally,  calculations  of  the  free  vibrations  of  complete  circular  plates 
with  clamped,  simply  supported,  and  free  boundaries  are  carried  out. 


2.0  THEORETICAL  FORMULATION 

The  stiffness,  mass,  and  load  matrices  for  the  finite  element  method  are 
most  easily  obtained  from  calculations  of  potential  and  kinetc  energies.  Once  the  dis¬ 
placement  functions  for  a  particular  finite  element  have  been  established,  they  may  be 
substituted  directly  into  the  potential  and  kinetic  energy  integrals. 

For  classical  bending  of  isotropic  plates,  the  potential  energy  for  an  element 
in  polar  co-ordinates  is  given  by  (page  346,  Ref.  3) 

82w  (l  9w  1  82w 

’’  3r2  lr  8r  r2  ao2 

2t  (1) 

t  if )  J  '  r  drd® 

where  q  is  the  external  loading  on  the  plate  element,  and  D  and  v  are  the  customary 
flexural  rigidity  and  Poisson’s  ratio,  respectively. 

The  kinetic  energy  for  an  element  in  polar  co-ordinates  is  simply 

T  ’  \  fo  fj*  <2> 

where  p  is  the  material  density  and  h  the  plate  thickness.  Making  the  usual  assump¬ 
tion  of  sinusoidal  dependence  on  time  yields 

2  rP  ,r  2 

T  =  f0  Jr*  Ph  w  r  dr  a*  (3) 

where  u>  is  the  frequency  of  oscillation. 

The  displacement  fuiicticr.s  for  the  two  finite  elements  developed  herein  are 
presented  in  Sections  2. 1  and  2.  2. 

2. 1  Circular  Sector  Finite  Element 


V-  =  f  fr 2  \d  f/sfw  1  3w  X  2-xf  .  211-1 

J0  \  i2  Lw  r  3r  r2  ao2  /  ‘ 


2(1  y)\r9r90 


The  first  finite  element  configuration  to  be  considered  is  shown  in  Figure  1. 


-  3  - 


K 

IS 


ESS 


The  element  is  a  sector  of  a  circular  plate  of  radius  r^  with  an  inv  luded  angle  $.  The 
corners  are  numbered  1  to  3,  as  shown. 


It  has  become  common  practice  in  developing  piate  bending  finite  elements  to 
use  at  least  the  displacement  and  two  slopes  as  generalized  co-ordinates  at  each  corner 
of  the  element.  This  allows  the  displacement  of  the  element  U>  vary  cubicallv  along  any 
edge.  Hence,  for  the  present  configuration,  this  suggests  the  use  of  at  least  the  follow¬ 
ing  nine  generalized  co-ordinates:  v.-rl9,  w-rl5,  w-j.  w  w  W0,  w  3,  w  3  and  W-  , 

where  w  =  Dw/5r,  and  w  ^  -  3w /dG.  The  single  number  subscripts  denote  the  comers 
of  the  element  shown  in  Figure  1.  and  w  is  the  slope  3w/9r  at  corner  1  in  the  direc¬ 
tion  of  corner  j,  for  j  =  2  or  3. 


Now  a  displacement  function  for  this  element  must  be  selected.  One  simple 
approach  to  developing  a  suitable  displacement  function  for  this  finite  element  configu¬ 
ration  appears  to  be  through  the  use  of  interpolation  formulas.  The  merit  of  this  ap¬ 
proach  for  developing  rectangular  finite  elements  has  been  clearly  demonstrated  in 
Inference  4.  i  he  interpolation  formulas  used  herein  are  some  of  the  Herir.ite  polyno¬ 
mials,  and  are 


H0i(x;  a)  =  *4j  (2x3-3ax2j-a3) 
a 

h02(x;  a)  =  ~  (3ax2-2x3) 
a" 

Hjj(x;  a)  =  ~  (x3-2ax2*a2x) 
a~ 

H,2(x;  a)  =  ~  (x3-ax2) 


Note  that  the  bracketed  superscript  1  that  usually  accomoanies  each  H..(.\)  has  been 

IJ 

omitted  for  simplicity.  Plots  of  these  functions  are  shovoi  in  Figure  2.  It  may  be  ob¬ 
served  that  the  Hn,  ai.d  Hno  functions  have  unit  value  at  x  =  0  and  x  =  a,  respec- 

tivelv;  zero  value  at  x  =  a  and  x  =  0,  respectively;  and  zero  slope  at  both  x  =  0 
and  x  -■  a.  On  the  other  hand,  the  Hj  j  and  H  functions  are  zero  at  both  x  =  0  and 

x  =  a,  but  have  unit  slope  at  x  =  0  and  x  =  a,  respectively.  Appropriate  combina¬ 
tions  of  these  tunctions  ?»-e  then  used  with  each  generalized  co-ordinate;  for  instance, 
for  a  rectangular  plate  (Ref.  4)  the  functions  that  accompany  a  typical  comer  displace¬ 
ment  are 


W!  H01(X;  a)  HQ1^  :  b) 


and  with  a  typical  slope  are 


w 


xl 


Hjjfx;  a)  H01(y:b) 


For  this  circular  sector  element  the  displacement  may  be  assumed  in 
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the  form 

w(r,  6)  =  wrl2  Hn(r;  rQ)  HQ1(0;  0)  +  Hn(r;  rQ)  11^(0;  0)  +  H01(r;  rQ) 

+  .  (5) 

It  may  be  seen  that  this  form  provides  the  proper  dependence  on  the  generalized  co¬ 
ordinates;  for  example,  the  first  term  provides  an  acceptable  dependence  on  w  *2  in 

that  H11(0;  rQ)  =  0,  d  H^/drfO;  rQ)  =  1  and  HQ1(0;  0)  =  1.  However,  upon  closer 

inspection,  it  is  seen  that  the  first  two  terms  of  expression  (5)  contain  linear  terms  in 
r,  which  represent  an  infinite  bending  moment  at  r  =  0.  This  therefore  precludes  the 
use  of  expression  (5). 

Since  there  is  no  way  to  avoid  linear  terms  in  r  and  still  retain  the  slope 
degrees  of  freedom  at  r  =  0,  different  interpolation  functions  must  be  found  for  the  0 
dependence.  It  is  logical  to  use  the  trigonometric  functions  sin  0  and  cos  0  in  these 
interpolation  functions,  since  r  sin  0  and  r  cos  0  are  harmonic  functions.  Hence,  the 
following  trigonometric  functions  are  found  to  be  satisfactory 

(pA9)  =  cos  0  -  cot  0  sin  0 
1  (6) 
<?2(0)  =  sin  0/sin  0 


Typical  plots  of  these  two  functions  are  shown  in  Figure  3  for  a  few  values  of  the  in¬ 
cluded  angle  0.  It  may  be  seen  that  <p  has  unit  value  at  0  =  0  and  is  zero  at  9  =  0, 

while  <?2  is  zero  at  0  =  0  and  has  unit  value  at  0  =  0.  Unfortunately,  the  slopes  are, 

in  general,  non-zero  at  either  0  =  0  or  0  =  0.  These  functions  have  the  additional 
disadvantage  that  their  distributions  between  0  =  0  and  0  are  not  independent  of  0. 
However,  it  is  very  difficult  to  find  any  better  functions  that  are  still  compatible  with 
the  linear  term  in  r.  Consequently,  the  functions  <p^  and  <?2  are  used  in  the  present 


derivation  by  assuming  that  w  12 


and  w 


rl3 


contribute  to  the  element  displacement  in 


the  following  way 


w(r,  0)  =  [wrl2  ^(0)  +  wrl3  <f>2<0)]  Hn(r;  rQ)  +  HQ1(r;  r0)  + 


(7) 


it  may  be  noted  that  the  form  of  equation  (7)  does  represent  well-behaved  bending  mo¬ 
ments  at  r  =  0. 


One  further  aspect  of  the  displacement  expression  (7)  should  be  noted  here. 
Differentiating  expression  (7)  with  respect  to  0  and  setting  0=0,  yields  the  circum¬ 
ferential  slope  along  the  element  edge  1-2.  It  may  be  seen  that  this  slope  depends  on 
w  since,  in  general,  3<po/80  is  non-zero  at  0  =  0.  This  means  that  the  normal 

I*  lt>  6 

slope  along  the  element  edge  1  -  2  is  not  independent  of  the  displacements  along  the 
element  edge  1-3.  Similarly,  it  may  be  seen  that  the  normal  slope  along  the  edge 
1  -  3  is  not  independent  of  the  edge  1-2.  This  implies  that  the  element  will  not  be  of 
the  conforming  type,  since  the  normal  slopes  will  not  be  continuous  along  the  radii  be¬ 
tween  adjacent  elements.  However,  it  must  be  emphasized  that  this  does  not  preclude 


^*a5fcS»£raS* 
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the  displacements  being  continuous  between  adjacent  elements. 

Finally,  since  the  element  cannot  be  conforming,  it  is  suffici  :nt  to  use  the 
nine  generalized  co-ordinates  suggested  earlier  and  to  assume  the  element  displace¬ 
ment  in  the  form 


+  w 


^2(0)]  Hll(r;  r0> 


w(r,  0)  =  wx  H01(r;  rQ)  +  [wfl2  ^(0) 

+  [Wr2  H01(0;  P)  f  Wr3  HO2(0;  0)]  H12(r;  r0) 
+  [W02  Hll^’  ^  +  W03  H12^0’  H02^r; 

+  [W2  H01<*;  «  +  W3  HO2<0;  «]  H02<r;  V 


(8) 


Note  that  w  is  a  cubic  function  of  r  throughout  the  element.  Furthermore,  along  the 
edge  1  -  2  (0  =  0),  this  cubic  variation  is  uniquely  defined  by  w^^,  wj»  wr2*  w2‘ 

The  analogous  result  also  holds  for  the  edge  1-3.  Hence,  the  displacement  will  be 
continuous  along  the  common  boundary  between  adjacent  elements,  provided  their  dis¬ 
placements  and  slopes  at  r  =  0  and  rQ  are  set  equal  on  that  boundary.  It  may  also  be 

seen  from  equation  (8)  that  along  the  circular  edge  2  -  3  (r  =  rQ),  w  is  a  cubic  function 

of  0  and  depends  only  on  wfl2,  w2,  w63,  and  wg. 

To  find  the  stiffness  matrix  for  e  element,  the  displacement  function  given 
by  equation  (8)  is  now  substituted  into  the  strain  energy  part  of  equation  (1),  and  the 
integration  over  the  element  is  carried  out.  Note  that  for  this  circular  sector  finite 
element,  the  integration  over  r  is  from  r^  =  0  to  r^  =  r^.  This  integration  yields 

the  strain  energy  in  the  form 


u  -  5  jxiiT  [Ki]  K 


(9) 


where  { >  is  the  column  vector 


jXl|  =  (wrl2’  wrl3*  wl-  wr2*  W«2-  w2-  wr3-  w03'  w3) 


(10) 


and  [kj]  is  the  9x9  stiffness  matrix  given  in  Table  I.  The  mass  matrix  is  found 
in  a  similar  manner  by  substituting  equation  (8)  into  equation  (3).  This  yields  the  kinetic 
energy  in  the  form 

t  -  4  hf  H  hi  <»> 

where  [Mj]  is  the  9  x  9  mass  matrix  given  in  Table  II.  The  consistent  load  matrix  for 
the  element  may  be  obtained  by  evaluating  the  wq  portion  in  the  integral  of  equation  (1). 
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2. 2  Angular  Sector  Finite  Element 

The  next  finite  element  configuration  to  be  considered  is  shown  in  Figure  4. 
The  element  is  a  sector  of  an  annular  plate  with  inside  radius  r^  outside  radius  r 2, 

and  an  included  angle  /?.  The  corners  are  numbered  1  to  4  as  shown. 

A  primary  objective  in  the  derivation  of  this  element  is  to  make  it  consistent 
and  compatible  with  the  circular  sector  element  of  the  previous  Section.  This  will  allow 
the  two  kinds  of  elements  to  be  put  together  for  the  solution  of  problems  involving  por¬ 
tions  of  circular  plates  that  include  the  centre  of  curvature.  Therefore,  it  will  be  suf¬ 
ficient  to  use  the  displacement  and  two  slopes  at  each  corner  as  the  generalized  co¬ 
ordinates  for  the  element.  This  will  allow  the  displacement  to  be  cubic  in  r  along 
edges  1-2  and  3-4,  and  cubic  in  0  along  the  curved  edges  1-4  and  2-3.  Hence, 
the  displacement  will  be  continuous  between  elements  with  a  common  boundary,  provided 
the  corner  displacements  are  made  equal  at  the  ends  of  that  boundary,  but  the  slope  nor¬ 
mal  to  that  boundary  will,  in  general,  not  be  continuous  across  it.  The  displacement 
vector  for  the  annular  sector  finite  element  then  becomes  the  twelve-term  column  vector 


where  the  subscripts  1  to  4  denote  the  corners  of  the  element  as  shown  in  Figure  4. 

The  displacement  function  for  this  element  is  assumed  in  the  form 
w(r,  8)  ~  aj  +  agr  +  ag0  +  a4r0  +  agr2  +  a ^8*  +  a?r20  +  agr02  +  agr3 

+  a1Q03  +  anr30  +  a12r03  (13) 

It  may  be  observed  that  this  function  is  similar  in  form  to  that  used  for  the  well-known 
twelve  degree  of  freedom  rectangular  plate  bending  element  (Ref.  1).  In  this  case,  the 
r  and  0  co-ordinates  just  replace  the  x  and  y  cartesian  co-ordinates,  respectively, 
used  for  the  rectangular  element. 

The  twelve  corner  displacements  used  in  equation  (12)  may  now  be  evaluated 
from  equation  (13).  Carrying  this  out  leads  to  the  matrix  relation 

K|  -  M  M  <i4> 

where  |a|  is  the  column  vector  of  polynomial  coefficients 


and  [t]  is  the  12  x  12  transformation  matrix  given  in  Table  III.  In  principle,  this 
transformation  matrix  may  be  inverted,  and  the  polynomial  coefficients  a.  may  be  de¬ 
termined  as  functions  of  the  comer  displacements  of  equation  (12).  The  a.  may  then 
be  substituted  into  equation  (13),  yielding  the  displacement  function  in  terms  of  the 
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corner  displacements.  However,  carrying  out  this  process  and  substituting  the  result¬ 
ing  displacement  function  into  the  energy  integrals  involves  an  amount  of  work  that  is 
impractical. 

An  alternative  approach,  which  requires  much  less  effort,  is  to  substitute 
equation  (13)  directly  into  the  energy  integrals.  Carrying  out  the  integrations  over  the 
element  yields  the  strain  energy  in  the  form 

u  - 1  iAf  M  lAl  <16> 

and  the  kinetic  energy  in  the  form 

T  *  4  !A!T  M  lAl  <l7> 

where  [k]  and  [m]  are  the  12  x  12  stiffness  and  mass  matrices  given  in  Tables  IV 
and  V,  respectively,  and  jA{  is  given  in  equation  (15).  Eliminating  )a|  from  equa¬ 
tions  (16)  and  (17)  by  means  of  equation  (14)  yields 

u  - 1  Kf  M  ta!  (18) 

and 

2 

T  -  !s2|t  Kl  lx2|  *i9> 

where 

T 

W  -  [t'1]  M  [t'1]  m 


and 

[m2]  -  [T-f  [m]  [t’1]  (21) 


are  the  required  12  x  12  stiffness  and  mass  matrices,  respectively,  for  the  element  in 
terms  of  the  comer  displacements.  The  calculations  involved  in  equations  (20)  and  (21) 
are  carried  out  on  the  computer  to  any  desired  accuracy.  The  consistent  load  matrix 
for  the  element  is  found  frcm  the  wq  integral  (see  eq.  (1)),  first  in  terms  of  the 
polynomial  coefficient  vector  )A|  and  then  transformed  for  the  corner  displacement 


vector  X9 1  by  premultiplying  by 


3.0  NUMERICAL  APPLICATIONS 

& 

The  two  finite  elements  developed  in  the  previous  Section  are  now  used  to 
solve  several  numerical  examples.  These  examples  serve  to  indicate  how  well  the 
elements  work  in  various  types  of  applications.  The  problems  considered  are  the 
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static<  nonsymmetric  deflections  of  an  annular  plate  and  a  complete  circular  plate,  both 
loaded  by  point  loads  and  the  free  vibrations  of  complete  circular  plates  with  simply 
supported,  clamped  and  free  boundaries. 

3. 1  Point  Load  on  Annular  Plate 

The  first  problem  considered  is  that  of  the  static,  nonsymmetric  deflection 
of  the  annular  plate  shown  in  Figure  5.  The  plate  is  clamped  along  the  inner  edge  (r  =  b) 
and  loaded  by  a  concentrated  force  at  the  outer  boundary  (r  =  a).  This  annular  plate 
configuration  is  chosen  because  it  may  be  modelled  with  only  the  annular  sector  finite 
elements,  and  hence  provides  a  check  of  those  elements  alone.  The  concentrated  load 
is  known  to  provide  one  of  the  severest  tests  of  the  finite  element  method.  The  problem 
is  also  a  convenient  one  to  use  here  because  it  has  an  exact  series  solution.  This  series 
solution,  which  is  given  on  pages  290  -  292  of  Reference  3,  was  programed,  and  the 
exact  displacement  was  calculated  for  comparison  with  the  finite  element  solutions. 

The  finite  element  calculations  were  carried  out  for  the  particular  case  dis¬ 
cussed  in  Reference  3:  a  =  1.  5,  b  =  1.  0,  and  v  =  0.  3.  Symmetry  about  the  dia¬ 
meter  containing  the  load  was  used,  so  that  only  one-half  of  the  plate  had  to  be  modelled. 
The  four  successive  refinements  of  finite  element  assemblages  shown  in  Figure  6  were 
used  in  the  calculations.  The  finite  element  results  for  the  displacement  under  the  load 
are  given  in  Table  VI,  along  with  the  number  of  degrees  of  freedom  used  in  each  calcu¬ 
lation.  The  exact  solution  is  also  presented  in  the  Table.  The  radial  and  circumferen¬ 
tial  displacement,  distributions  obtained  from  the  4  x  24  grid  of  finite  elements  are  plot¬ 
ted  in  Figure  7  along  with  the  exact  ones. 

Table  VI  shows  that  the  displacement  under  the  load  does  not  converge  mono- 
tonically  with  refinement  in  finite  element  size.  This  is  not  surprising  because  the 
finite  elements  are  not  of  the  conforming  type.  However,  the  displacements  shown  in 
Table  VI  are  still,  quite  accurate  and  appear  to  be  converging  towards  the  exact  value 
as  the  element  size  becomes  quite  small.  It  should  be  noted  that  even  though  the  dis¬ 
placement  under  the  load  predicted  by  the  1x6  grid  of  elements  was  quite  close  to  the 
exact  value,  the  circumferential  displacement  distribution  away  from  the  load  was  very 
inaccurate.  This  is  shown  in  Table  VII,  where  numerical  results  are  given  for  the  dis¬ 
placement  of  the  free  edge  at  various  circumferential  locations.  These  results  serve 
to  illustrate  the  type  of  convergence  obtained  away  from  the  point  load.  At  0  =  180°, 
this  convergence  appears  to  be  quite  slow,  but  it  may  be  noted  that  the  displacement 
there  is  six  orders  of  magnitude  smaller  than  under  the  load  ( 0  =  0).  In  summary,  it 
appears  that  the  annular  sector  finite  elements  work  very  well  for  the  type  of  problem 
considered  here. 

3. 2  Point  Load  on  Circular  Plate 

The  next  problem  considered  is  that  of  the  static,  nonsymmetric  deflection  of 
the  complete  circular  plate  shown  in  Figure  8.  The  plate  is  clamped  around  the  outer 
boundary  (r  =  a)  and  is  loaded  by  a  single  concentrated  force  P  at  r  =  b.  The  exact 
solution  for  this  problem  is  given  by  equation  (197)  on  Page  293  of  Reference  3. 

The  modelling  of  this  plate  requires  the  use  of  both  the  annular  and  circular 
sector  finite  elements  developed  herein.  Symmetry  about  the  diameter  containing  the 
load  is  again  used,  so  that  only  one-half  of  the  plate  is  solved.  The  numerical  calcula¬ 
tions  are  carried  nut  for  a  =  1.  0,  b  -  0.5,  and  v  =  0.3,  using  the  three  finite 
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element  assemblages  shown  in  Figure  9. 


The  numerical  results  for  the  displacement  under  the  load  are  presented  in 
Table  VIII,  along  with  the  numbers  of  degrees  of  freedom  used  in  each  calculation.  The 
displacements  predicted  for  a  few  other  positions  on  the  plate  are  given  in  Table  IX. 
Radial  and  circumferential  displacement  distributions  obtained  from  the  6x8  grid  cal¬ 
culations  are  shown  in  Figure  10.  All  these  results  exhibit  a  high  accuracy,  but  unfor¬ 
tunately  the  convergence  is  somewhat  disappointing.  For  example,  the  2x4  and  4x6 
grid  predictions  for  the  displacement  under  the  load  are  only  in  error  by  about  one  per¬ 
cent,  and  appear  to  be  converging  in  that  the  4x6  grid  answer  is  more  accurate  than 
the  2  x  4  one.  However,  the  6  x  8  grid  prediction  is  in  error  by  about  two  percent,  and 
the  sign  of  this  error  is  different  from  that  of  the  former  ones.  Since  the  results  in  the 
previous  Section,  using  the  annular  sector  finite  elements  by  themselves,  appeared  to 
be  well-behaved,  it  is  suspected  that  the  erratic  convergence  behaviour  observed  here 
is  probably  caused  by  the  circular  sector  finite  elements.  As  pointed  out  in  Section  2. 1, 
the  trigonometric  interpolation  functions  used  in  developing  these  elements  changed 
shape  as  the  included  angle  ft  was  changed.  However,  in  the  limit  of  very  small  in¬ 
cluded  angle,  these  functions  approach  a  constant  straight  line  shape.  Hence,  it  is  ex¬ 
pected  that  this  erratic  convergence  behaviour  would  disappear  if  the  finite  element  as¬ 
semblages  were  increasingly  refined.  In  summary  then,  even  though  the  convergence 
is  somewhat  erratic,  the  numerical  accuracy  is  still  quite  satisfactory  for  most  appli¬ 
cations. 

3.  3  Vibration  of  Circular  Plates 

The  free  vibrations  of  complete  circular  plates  with  clamped,  simply  sup¬ 
ported,  and  free  outer  edges  are  now  analyzed  with  the  present  finite  elements.  These 
final  applications  serve  to  check  out  the  consistent  mass  matrices  for  the  elements. 
The  calculations  are  carried  out  for  a  plate  radius  of  unity  and  Poisson’s  ratio  equal  to 
0.33,  using  the  three  finite  element  assemblages  shown  in  Figure  11.  As  indicated  in 
the  Figure,  the  3  x  12  grid  solution  was  separated  into  two  sub-problems  by  using 
symmetry  about  a  diameter.  This  was  done  in  order  to  keep  the  eigenvalue  problems 
to  a  manageable  size. 


The  non-dimensional  frequency  predictions  for  the  first  few  vibration  modes 
are  presented  ir.  Tables  X,  XI,  and  XII  for  the  three  problems.  The  exact  solutions 
are  taken  from  Reference  5.  It  may  be  seen  that,  in  some  cases,  two  different  eigen¬ 
values  were  obtained  for  apparently  the  same  vibration  mode.  Upon  close  inspection 
of  the  predicted  mode  shapes  it  was  found  that  the  predicted  frequencies  depended  on 
whether  or  not  the  radial  nodal  lines  fell  between  or  along  the  finite  element  junctions. 
For  example,  in  the  2x8  grid  predictions  of  Modes  Seven,  for  the  clamped  and  simply 
supported  plates,  and  Mode  Five  for  the  free  plate,  the  nodal  lines  were  along  the  ele¬ 
ment  junctions  for  the  low-  eigenvalues  and  mid-way  between  them  for  the  high  ones. 
The  same  trend  was  found  to  hold  for  the  other  non-repeated  eigenvalues.  This  peculiar 
effect  may  again  be  associated  with  the  circular  sector  finite  elements  and  their  special 
interpolation  functions.  It  may  be  noted  that  the  differences  between  the  non-repeating 
eigenvalues  appear  to  decrease  as  the  modelling  is  refined,  until  the  final  predictions 
are  quite  close  to  the  exact  values.  Again,  very  little  can  be  said  about  the  particular 
convergence  of  these  calculations  because  the  elements  are  not  of  the  conforming  type. 
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4.0  CONCLUDING  REMARKS 

A  finite  element  metnod  for  plate  bending  elements  in  polar  co-ordinates  has 
been  presented.  Two  elements  were  developed  -  one  in  the  form  of  a  circular  plate 
sector  and  the  other  in  the  form  cf  an  annular  plate  sector.  The  application  of  these 
elements  to  several  static  and  dynamic  numerical  examples  indicated  the  accuracy  to 
be  expected  with  the  method. 

These  elements  may  be  combined  with  rectangular  elements  for  the  analysis 
of  a  variety  of  plate  bending  problems  heretofore  considered  unsolvable. 
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TABLE  I 


STIFFNESS  MATRIX  JkJ/D  ror  CIRCULAR  SECTOR  FINITE  ELEMENT 
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TABLE  HI 

TRANSFORMATION  MATRIX  [t]  BETWEEN  CORNER  DISPLACEMENTS  AND 
POLYNOMIAL  COEFFICIENTS  FOR  ANNULAR  SECTOR  FINITE  ELEMENT 
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DISPLACEMENT  UNDER  CONCENTRATED  LOAD  FOR 


TABLE  Vin 


DISPLACEMENT  UNDER  CONCENTRATED  LOAD  FOR 
CLAMPED  CIRCULAR  PLATE 


FINITE 

ELEMENT 

GRIDS 


NUMBER  OF 
DEGREES  OF 
FREEDOM 


DISPLACEMENT 
UNDER  LOAD  P, 
wD/P 


2X4 

4X6 

6X8 


18 

63 
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0. 0113155 
0. 0112715 
0. 0109738 


EXACT  SOLUTION 
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TABLE  EX 

VARIOUS  DISPLACEMENTS  wD/P  x  103  OF 
CLAMPED  CIRCULAR  PLATE 
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TABLE  X 

NON-DIMENSIONAL  FREQUENCIES  y'ph/D  wR2  FOR 
CLAMPED  CIRCULAR  PLATE,  R  =  1.0 
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MODES 

FINITE  ELEMENT  GRIDS 
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MODES 


FINITE  ELEMENT  GRIDS 


3  x  12 


EXACT 

SOLUTIONS 


5. 3349 
9.  8543 


5.  9823 


5. 9071 
5. 9405 


■■CO 


9.  0941  9.  0818  8.  9795 


A7^ 


ran 


25. 155  13.  519  12.  982 


12.  25 


N 

5 

H 

20.  916  20. 160  20.  243 


23. 597 
25.  565 


39.  919 


23. 016 


34. 184 
34. 442 


33. 269  I  36. 014 


-  CO 


38.721  38.  116 


0 


x 


o 


FIG.  2  -  POLYNOMIAL  INTERPOLATION  FUNCTIONS 
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ANGULAR  CO-ORDINATE  6 


FIG.  3:  TRIGONOMETRIC  INTERPOLATION  FUNCTIONS 
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FIG. 5  POINT  LOADED  ANNULAR  PLATE  CONFIGURATION 
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RAOIAl  DISTRIBUTION 


FIG. 7-  DISPLACEMENT  DISTRIBUTIONS  FOR  CLAMPED-FREE 
ANNULAR  PLATE  UNDER  POINT  LOAD 


4  x  6  GRID 


P/2 

6x8  GRID 


FIG. 9:  FINITE  ELEMENT  ASSEMBLAGES  FOR 
ONE-HALF  CIRCULAR  PLATE 
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FIG. II :  FINITE  ELEMENT  ASSEMBLAGES  FOR 
CIRCULAR  PLATE  VIBRATIONS 


